Elongation and fluctuations of semi-flexible polymers in a nematic solvent 
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We directly visualize single polymers with persistence lengths ranging from ip = 0.05 to 16 fim, 
dissolved in the nematic phase of rod-like fd virus. Polymers with sufficiently large persistence 
length undergo a coil-rod transition at the isotropic-nematic transition of the background solvent. 
We quantitatively analyze the transverse fluctuations of semi-flexible polymers and show that at 
long wavelengths they are driven by the fluctuating nematic background. We extract both the 
Odijk deflection length and the elastic constant of the background nematic phase from the data. 

PACS numbers: 61.30.-v, 64.70.Md, 82.35.Pq 



Polymer coils in solution exhibit a variety of confor- 
mational and dynamical behaviors depending on many 
factors, including polymer concentration, polymer stiff- 
ness, solvent quality, solvent flow, and mechanical stress. 
Exciting recent experiments in this field have focused on 
disentanglement of single biopolymers in isotropic solu- 
tions as a result of applied forces and solvent flow 1] and 
on transport of single biopolymers through networks of 
barriers 0] , which is critically affected by conformational 
dynamics of the polymers. In this Letter, we explore con- 
formations of polymer coils in anisotropic solutions. In 
particular, we present the first direct experimental obser- 
vations of isolated semi-flexible polymers dissolved in a 
background nematic phase composed of aligned rod-like 
macromolecules. We show by direct visualization that 
semi-flexible biopolymers dissolved in the nematic phase 
assume an elongated rod-like configuration aligned with 
the background nematic director. The coil-rod transition 
of the biopolymer can thus be induced by causing the sol- 
vent to undergo an isotropic-to-nematic (I-N) transition 
by increasing the concentration of its constituent rods. 
We quantitatively explore the fluctuations of these semi- 
flexible polymers and find they cannot be described by 
a theory which treats the nematic background as a fixed 
external field 

Mixtures of semi-flexible polymers in lyotropic nematic 
suspensions exemplify an emerging class of complex flu- 
ids - hyper-complex fluids, for example nematic elas- 
tomers and nematic emulsions wherein two or 
more distinct components are combined to create sys- 
tems that exhibit novel physical properties and func- 
tions. Understanding the polymer- nematic system may 
lead to new ideas about how to achieve high alignment of 
biopolymers that is complementary to existing methods 
of DNA alignment j^j. Furthermore, since many biopoly- 
mers such as the actin filaments within the sarcomere and 
neurofilaments within the axon reside in an anisotropic, 
nematic-like environment (7J , our investigation may shed 
light on organization mechanisms within the cell. 

We have used fluorescence microscopy to study four 
different biopolymers in isotropic and nematic colloidal 



TABLE I: The contour length L, the persistence length £j 
and the diameter a of bio-polymers in our experiments. 



Polymer 


L [jJ.m] 


£p [fim] 


a [nm] 


Ref. 


A-DNA 


16 


0.05 


2 


[18] 


neurofilament 


2-10 


0.2 


10 


\m 


wormlike micelles 


5-50 


0.5 


14 




F- act in 


2-20 


16 


7 


[20] 


fd virus 


0.9 


2.2 


7 


[13j 



suspensions. This approach yields new information about 
dynamics and defects not readily accessible to traditional 
probes such as x-ray or neutron scattering In addi- 
tion, we have developed a rotationally-invariant free en- 
ergy for a single semiflexible polymer in a nematic matrix 
which generalizes the work in [3, E| , and enables us to 
extract the Odijk length and the elastic constant of 
the liquid crystal. These first direct measurements of 
the Odijk deflection length, A, allow us to quantify the 
length scale over which the polymer wanders before it is 
deflected back by the nematic director. 

Our experiments employ an aqueous solution of rod- 
like fd viruses as a background nematic li quid cry stal. 
This system has been studied extensively Ed; 13 ' and 
its phase behavior is well described by the Onsager theory 
for rods with hard core repulsion |l5j | . Another advantage 
of this system is its compatibility with most biopolymers. 
We use four different semi- flexible polymers, whose physi- 
cal parameters are listed in Table To directly visualize 
the polymers dissolved in the nematic background, we 
fluorescently labelled each polymer: DNA was labelled 
with YOYO-1 (Molecular Probes, Eugene OR), neurofll- 
aments with succinimidyl rhodamine B , F-actin flla- 
ments with rhodamine-phalloidin (Sigma, St. Louis MO), 
and wormlike micelles with PKII26 dye (Sigma, St. Louis 
MO) which preferentially partitions into the hydropho- 
bic core of the micelle. Since DNA, neurofilaments, and 
actin are all negatively charged, we expect that they are 
stable in a suspension of negatively charged fd viruses. 
Wormlike micelles are sterically stabilized with a neutral 
PEO brush layer, which does not interact with fd virus 
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or other proteins 0, 0| . 

Bacteriophage fd was grown and dialyzed against a 
phosphate buffer (150 mM KCl, 20 mM phosphate, 2 
niM MgCl2, pH=7.0) (l^l- Samples were prepared by 
mixing a small amount of polymer with fd solution at 
different concentrations and were placed between a cov- 
erslip and a glass slide. A chamber with a thickness of 
~ 50 /im was made by using a stretched parafilm as a 
spacer. Samples sealed with optical glue (Norland Prod- 
ucts, Cranbury, NJ) were allowed to equilibrate until no 
drift was visually detectable. To reduce photoblcach- 
ing, we added anti-oxygen solution (2 mg/ml glucose. 
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FIG. 1: Images of fluorescently labelled biopolymers in the 
isotropic (left) and nematic (right) phase of /rf virus. Figures 
(a)-(d) are, respectively, the images of actin, wormlike mi- 
celles, neurofilaments, and DNA. The polymers in an isotropic 
solution are confined by a thin chamber thus making the sam- 
ples quasi two dimensional, (e) A sequence of images illus- 
trating an actin filament escaping from a hairpin defect. The 
scale bar is 5 pm. (f) A schematic of a biopolymer in the 
background nematic field; the conformation of the polymer 
is parameterized by R(z) = {rx(z\ry{z),z) . The nematic 
director points along the z-axis. 



360 U/ml catalase, 0.25 vol% mercaptoethanol, 8 U/ml 
glucose oxidase). All samples were imaged with a fluo- 
rescence microscope (Leica IRBE) equipped with a lOOx 
oil-immersion objective and a 100 W mercury lamp. Im- 
ages were taken with a cooled CCD camera (CoolSnap 
HQ, Roper Scientific), which was focused at least 5 ^m 
away from the surface to minimize possible wall effects. 

Figure ^ displays a series of pictures that summarize 
our qualitative observations. In the nematic phase of 
/d, F-actin (Fig. ^), wormlike micelles (Fig. ^), and 
neurofilaments (Fig. ^) are highly elongated, having a 
rod-like shape. By contrast, the same filaments dissolved 
in an isotropic phase crumple into more compact random 
coils. Just above the I-N transition, actin filaments and 
worm-like micelles form hairpin defects These hair- 
pins exhibit interesting dynamics as shown in Fig.^Ji f^nd 
will be explored by us in detail elsewhere. DNA dissolved 
in /d nematic behaves qualitatively differently (Fig. ^); 
it forms a slightly anisotropic droplet. Each droplet con- 
tains many DNA molecules and, with time, these droplets 
coalesce into a larger droplet. Thus, even at a very 
low concentration, DNA separates from the /rf nematic. 
Taken together, these observations suggest that the per- 
sistence length of the polymer is important in determin- 
ing its solubility in the nematic liquid crystals: DNA has 
a small 1^ and is insoluble, unlike the other stiffer poly- 
mers in our experiments. This insolubility may be re- 
lated to the entropy-driven phase-separation of a system 
of bidisperse rigid rods if their lengths and/or diameters 
are sufficiently dissimilar j21j. This theory, however, has 
not been extended to the case of semi-flexible polymers. 

The large contour lengths of actin filaments and worm- 
like micelles make them suitable for further quantitative 
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FIG. 2: The order parameter of actin filaments iSacu-n) vs. the 
order parameter of the background /d nematic (Sfd)- Dashed 
line is a guide to the eye. The contour length of actin fila- 
ments is 15 /im or higher. The values of Sfd at different fd 
concentration have been measured in Ref. p3 |. Inset: The 
orientational distribution function (ODF) of actin filaments. 
The ODF is well approximated by a Gaussian distribution for 
a wide range of concentrations. 
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analysis. We focus on the fluctuations of filaments in a 
background nematic that is free of both defects and dis- 
tortions. A series of 50 to 100 images where taken with a 
few seconds between each image to ensure that statisti- 
cally independent configurations were sampled. For each 
fd concentration, ten filaments were analyzed. The con- 
formation of each polymer was reconstructed by manu- 
ally marking the end points. (Note that we parameterize 
the transverse deviations of the polymer from the z axis 
by the 2-component vector r{z), as shown in Fig. ^.) 
An intensity profile along the x direction for each value 
of z was extracted. By fitting this intensity to a Gaus- 
sian, we obtained sub-pixel accuracy for rx{z). We first 
extracted the orientational distribution function (ODF) 
from our data. Since our images are two dimensional pro- 
jections of the polymer fluctuating in three dimensions, 
the x-component of the tangent vector is measured by 
tx{z) — drx{z)/dz K, dx(z). The ODF is obtained by 
creating a histogram of 9^ at different positions along the 
contour length for a time sequence of 50-100 images. A 
typical ODF is plotted in Fig. 12^-; it is well approximated 
by a Gaussian distribution. 

Next, we compute the order parameter of the poly- 
mer deflned by: S = dz{3{t{z) ■ zf - 1)/{2L) = 
1 — 3 ( (0) ) , where in the last equality, we have used 
{t1.{z)) — {ty{z)). In Fig. 121 we plot S for actin as a 
function of the background nematic order parameter. It 
is interesting to observe that SacUn is signiflcantly higher 
than Sfd- In order to check that the difference in the 
alignment between actin and fd molecules is due to their 
different contour lengths, we measured SacUn for different 
contour lengths of actin filaments, as shown in Fig.O As 
the contour length of actin decreases, SacUn approaches 
Sfd, as expected intuitively. These observations are qual- 
itatively consistent with the Onsager theory for a bidis- 
perse mixture of rod-like particles with different lengths 
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FIG. 3: SacUn vs. contour lengths of actin. The concentra- 
tions of the background nematic fd are 41 mg/ml (squares, 
Sfd = 0.75) and 28 mg/ml (circles, Sfd = 0.855). Dashed 
Unes are a guide to the eye. 



considered in Ref. |22|. This theory predicts the order 
parameter of long rods will be higher than the order pa- 
rameter of the background nematic of shorter rods. 

Finally, we measured the tangent-tangent correlation 
fimction (TTCF) {tx{z') tx{z' + z) ) for wormhke micelles 
dissolved in fd virus with concentration 40 mg/ml and 
above (Fig. At low fd concentrations, the fluctua- 
tions of worms are large as evidenced by visual observa- 
tion of spontaneous formation and dissolution of hairpin 
defects. In this regime, the measured TTCF does not de- 
cay uniformly. We thus focus our analysis on the regime 
where the background order parameter is very high, and 
the amplitude of the polymer fluctuations is small. This 
makes our data suitable for comparison with the theoret- 
ical model outlined below. 

The fluctuations of a semi-flexible polymer in a ne- 
matic phase may be described by the free energy [Toll23| : 



F = 



knT 







j dzl^tj, 






dz 



+ r[ti(z)-5n(0,z)]' 



i j d^xK{V6n) 



where kg is the Boltzmann constant, T is the tempera- 
ture, £p is the persistence length of the semi-flexible poly- 
mer, r is the strength of the coupling of the polymer to 
the background nematic fleld, Sn is the local direction of 
the fluctuating nematic fleld, and K is the nematic elas- 
tic constant. Note that and 6n are two-dimensional 
vectors in the plane perpendicular to the average direc- 
tor. It is straightforward to compute ( tx{z') tx{z' + z) ) 
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FIG. 4: The x component of the tangent-tangent correlation 
function for wormlike micelles measured at three different fd 
concentrations (cfd)- With increasing fd concentration the 
magnitude of the correlation function decreases. The solid 
lines are theoretical curves generated from Eq. with the 
best-fit parameters listed in Table ITTl Inset: TTCF for the 
lowest concentration of the fd virus. The dashed and dotted 
lines are respectively the contributions of the first and second 
term in Eq. (|5J. The data points below 0.5 /im are unreliable 
and have been excluded from the fitting ,24,1. 
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TABLE II: The Odijk deflection length A, the elastic con- 
stant of the background nematic K, and the coupling con- 
stant between wormlike micelles and background nematic T 
for different fd concentrations obtained from the fits shown in 
Fig-E] The best-fit value of (.p of wormlike micelles is 1.5 ^m. 



Cfd [mg/ml] 


A [nm] 


K [10-" dyne] 


r [l/fim] 


39 


0.18 


1.9 


46 


51 


0.13 


2.4 


88 


97 


0.06 


2.8 


416 



from Eq. 



(2) 



X / dx 



cos(xz/A) log(l -I- L'Va;^) 



(l+a;2)[l + ^2+ r3|l0g(l+^2/^2)] 



where D — 27rA/a is related to the molecular cutoff which 
we assume to be the diameter of the polymer (a ~ 10 
nm). The first term in Eq. jSJ describes the fluctua- 
tions of a semi-flexible polymer in a static external field, 
with a decaying length set by A = y^ip/T. The second 
term describes the fluctuations of the polymers driven by 
the tight coupling to the fluctuations of the background 
nematic field. Note that this term generalizes that of 
Ref. 10], in that it includes the back reaction of the 
stiff polymer on the nematic fluctuations. Since it de- 
cays approximately as a power law, we expect at large 
lengthscales the fluctuations of the polymer are always 
dominated by the nematic fluctuations. 



Figure0]shows our measured TTCF along with the flt- 
ted curve of Eq. (jSJ. Overa ll, g ood agreement is obtained 
at distances above 0.5 fxm |2J|. At these distances, most 
of the fluctuations of the worms are driven by the tight 
coupling to the background nematic fleld, coming from 
the second term in Eq. ((2Jl. The best-flt value of £p is 
found to be 1.5 ^m, somewhat higher than that obtained 
in previous measurements 17]. From the fits to the data, 
we extract the values of the Odijk deflection length A, K, 
and r, as listed in Table UTI We observe that with increas- 
ing /rf concentration, A decreases, while K and F increase, 
as one would intuitively expect. Finally, we note that 
the values for K are in agreement with previous mea- 
surements of twist elastic constant K22 = 3 • 10~® dyne 
for fd samples prepared under similar conditions |25| . 

In conclusion, we have shown that semi-flexible poly- 
mers with large enough persistence lengths assume a rod- 
like conformation when dissolved in a nematic solvent. 
Using image analysis, a full nematic orientational dis- 
tribution function was measured. In addition, we have 
shown that fluctuations of the polymer are driven pri- 
marily by the fluctuations of the background nematic 
fleld. Direct visualization of individual polymer yields 
valuable new information about the behavior of polymer 
chains in anisotropic solvents. 
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